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Abstract 

We consider the interactions between two uncharged planar macroscopic surfaces immersed 
in an electrolyte solution which are induced by interfacial selectivity. These forces are taken into 
account by introducing a depletion free-energy density functional, in addition to the usual mean- 
field Poisson-Boltzmann functional. The minimization of the total free-energy functional yields 
the density profiles of the microions and the electrostatic potential. The disjoining pressure is 
obtained by differentiation of the total free energy with respect to the separation of the surfaces, 
holding the range and strength of the depletion forces constant. We find that the induced 
interaction between the two surfaces is always repulsive for sufficiently large separations, and 
becomes attractive at shorter separations. The nature of the induced interactions changes from 
attractive to repulsive at a distance corresponding to the range of the depletion forces. 

1 Introduction 

Electrostatic interactions often play an important role in a variety of different systems, ranging 
from biological membranes to chemical industrial paint ingredients. In some cases, it provides 
the underlying mechanism for the stabilization of mesoscopic systems against flocculation and 
precipitation. When two macroscopic charged surfaces approach one another, the result is usually 
a repulsive force, which inhibits a further approach. For two flat charged plates, this effect can be 
understood in a physical picture in terms of the osmotic pressure generated by the difference of the 
ion concentration in the region between the two approaching surfaces and the electrolyte-reservoir 
concentration. On the other hand, attractive interactions, which lead to aggregation or fusion, are 
sometimes a desirable feature. This is the case, for example, in the adhesion and fusion of vesicles 
and membranes or in environmental sewage treatment. Furthermore, some experiments [jl], ^, |3|, [3], [|] 
and simulations || 0, § indicate that, for small separations and high surface-charge densities, 
two like-charged polyions can indeed attract. 

From the theoretical point of view, several distinct mechanisms leading to attractive interactions 



have been proposed, which are based on charge fluctuations [10, 11, 12], strong positional charge 



correlations [13, |I4], [l5| ], anisotropic hypernetted chain calculations [16] or strong bulk-counterions 
correlations [17, pLq] . Very recently a unified treatment taking into account quantum fluctuations and 
structural correlations of the Wigner crystals formed by the condensed counterions onto the charged 
surfaces has been proposed [O]. Although the bare Coulomb force between two macroscopic sur- 



faces is always repulsive, correlations and/or fluctuations can induce attractive interactions, which 
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occasionally may overcome the electrostatic repulsion between the two equally charged surfaces. 
Correlations, which are entirely neglected within the mean-field Poisson-Boltzmann (PB) approxi- 
mation (Gouy-Chapman theory^, ^]), are believed to be essential ingredients for the appearance 
of attractive interactions. Thus, most proposed mechanisms which lead to attraction always include 
a non-mean-field effect. 

In this work we propose a new mechanism for attraction between two identical plates. In 
contrast to the previous theoretical pictures, this mechanism is entirely at the mean-field level. 
However, non-pure electrostatic forces are taken into account by including depletion forces — for 
example, those associated with finite ionic radii — acting on one of the ion species surrounding 
the plates. For simplicity, we consider the case in which the identical plates are uncharged and 
infinitely large. By considering uncharged plates we can discern the effect of the depletion forces 
separately from the usual electrostatic mean-field repulsion, which indeed turns out to be entropic 
and not strictly electrostatic. If we treat the surface- charged case, we are not able to separate the 
two contributions. Due to the simplicity of the model, it is possible to derive explicit, analytical 
expressions for all thermodynamical properties, including the disjoining pressure. 

The remainder of this paper is organized as follows. In Section 2 the model is introduced and 
the general equations are obtained. Section 3 is devoted to solving the generalized PB equations for 
the non-overlapping regime, when the depletion zones associated with the two plates do not overlap. 
The solution to the generalized PB equations for the overlapping regime, when the depletion zones 
associated with the two plates do overlap, is obtained in Section 4. Some concluding remarks are 
presented in Section 5. The closed analytical expression for the disjoining pressure is obtained in 
the Appendix. 



2 Definition of the model 

We shall consider two uncharged macroscopic surfaces immersed in a symmetric 1:1 electrolyte 
within mean-field theory. The system is modeled by two planar, infinitely thin rigid and uncharged 
surfaces, separated by a distance h, in contact with a monovalent salt reservoir of bulk concentration 
uq. A Cartesian coordinate system is chosen so that the surfaces are located at the x = ±h/2 planes, 
in such a way that the x axis is perpendicular to the surfaces. At the mean-field level the microions 
are treated as an inhomogeneous ideal gas, with local number densities n+{x) and n_(cc) for the 
positive and negative ions, respectively. We assume that, due to some non- electro static depletion 
mechanism, these local densities become inhomogeneous in the region close to the infinite plates. 
These inhomogeneities are governed by the reduced (total) free-energy functional (per unit area), 
f = /3f, where P = l/k B T, 

/ = /depletion + fpB , (1) 

which we split into two terms. The first term of (|l|) corresponds to a non- electrostatic depletion 
free energy (per unit area), 

oo 

/depletion = e / &Xn + {x) 
— oo 

where e is a depletion-strength parameter and w s (C) can be considered a normalized external (non- 
electrostatic) potential with a finite short range, s. This term breaks the original degeneracy 
between cations and anions, penalizing positive particles that are closest from a distance s to the 



w. 



x + 



+ w s 



(2) 
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surfaces. It mimics, for example, the effect of different sizes for the microions. Smaller negative 
ions are allowed to come in direct contact with the neutral surfaces, whereas the positive particles, 
by their larger size, are held apart from an effective distance s, related to their sizes. The effect 
of this term on the system is to yield an excess of anions in the region surrounding the plates, 
leading to an inhomogeneity of the local densities of microions in the vicinity of the uncharged 
plates. Thus, although the surfaces are themselves neutral, this imbalance of microions gives rise 
to a non-vanishing electric field. To allow analytical calculations, we shall hereafter assume that 
w s has the step-function form, 

' 0, |5| > a, 

In the limit s — > 0, the function w s (£) corresponds to the Dirac delta function, 5(£) = limu; s (£). 

The second term of @), fpB, represents the reduced bulk excess PB free-energy functional (per 
unit area), 



(3) 



PB 



dx jri-f (z) (in A 3 n + (x) 



1 )+n. 



(x) (in A 3 n_(x) - lj 



+ 



(x) [n+(x) - n_(x)] - Pn [n + (x) + n_(x)] + f3U 



dx < n + {x) In [n+(x)/no] + n-(x) In [n_(x)/no] + -<i>{x) [n + (x) — n_(x)] 



- [n+(x) - 2n ] j , (4) 

where A is an arbitrary length scale. The electrochemical potential and the reference pressure 
were set, respectively, to (3/j, = In (A 3 no) and ftHo = 2no, since the system is in electrochemical 
equilibrium with the infinite salt reservoir. The reduced electrostatic potential, <f>(x) = (3eip{x), 
where e is the proton charge and ip(x) is the electrostatic potential, satisfies the Poisson equation, 



d 2 (f>(x) 
dx 2 



-4:Tr£ [n+(x) — 7i— (x)] , 



(5) 



where I = (5e 2 /D is the Bjerrum length and the solvent is treated as a continuum of dielectric 
constant D. 

Minimization of the reduced total free-energy functional f[n+(x),n-(x)] with respect to the 
number densities, 



5f[n + (x),n_(x)] 

8n + (x) 
6f[n+(x),n-(x)] 



In [n+(x)/no] + <j>(x) + e 
In [n_(x)/no] — 4>{x) = 



h 



x + - ] + w s \x — — 



h 



0, 



8ri-(x) 

leads to the Boltzmann distribution for the optimum microion profiles, 



n+(x) = no exp 



-(f)(x) 



h 

Ul\ j X + - 



< ic, ! x - — 



7i-(x) = noexp[0(x)] 



(6) 
(7) 

(8) 
0) 
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Replacing (|8|) and @ into the Poisson equation (||) leads to a generalized PB equation, 



d 2 d>(x) k 2 



dx 2 



y {exp [<j>{x)] - exp 



h\ ( h 



-4>{x) - ew s ( x + - j - ew s I x - - 



(10) 



where k = ^/&KnQi is the inverse of the Debye screening length. 

The appropriate boundary conditions are the vanishing of the electrostatic potential and the 
electric field at infinity, 

4>(x ±00) = <t>'(x ±00) = 0; (11) 
the vanishing of the electric field at the midplane (x = 0), 

(j)'( x = 0) = 0; (12) 

the continuity of the electrostatic potential and the electric field across the planes located at x = 
±f±s, 

( j ) (x]±^±s) = ^xi±^±s), (13) 
<P'(x]±^±s) = <p'(xl±^±s), (14) 

where (f)(x f y) = lim 4>{x) and <p(x I y) = lim <fi(x). The boundary conditions (]l3| ) and (ll^ ) are 

based on the fact that the charge distribution, which appears on the right-hand side of the Poisson 
equation contains just a finite jump at the planes x = ±| ± s. 

By symmetry we have (j)( x ) = <K — x ) an d we need only to consider the positive x axis. Because 
of the non-electrostatic depletion of cations around the surfaces located at x = ±|, the electrostatic 
potential (j){x) is always negative, since there is an effective excess of anions around the surfaces. 
We shall consider two regimes separately, namely, the non-overlapping regime {h > 2s) and the 
overlapping regime (h < 2s). 



3 The non-overlapping regime, h > 2s 

In the non-overlapping regime, which occurs when the separation between the surfaces is larger 
than the range of the depletion forces, h > 2s, the depletion zones associated with the two interfaces 
do not overlap, and the generalized PB equation reads 



k 2 sinh d>(x), for < x < ft — s and x > 4 + s, 

— — N - - 2 - 2 ( 15 ) 

dx z y e ~ 2a K 2 sinh [<f>{x) + 2q] , for § - s < x < § + s, 

where we introduced the parameter a = e/8s. 

Using the identity d ^2 - = 5 d [fJ ; the non-linear second-order differential equation represented 



by (|lq ) can be analytically integrated. Introducing the midplane electrostatic potential, 4> m = 
4>{x = 0), and the internal and external electrostatic potentials in the vicinity of the interface at 
x = |, (f) < = 4>{x = I — s) and 0> = cj){x = | + s), the solution which satisfy the boundary 
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conditions ([TT]) and @ can be written explicitly as 

kA [<j) m ,(fr(x)} , 



for < x < 4 - s, 



</>'(x) 



Ksign (x — xi) e a A Q [0i, </>(x)] , for 



-2k sinh 



#■0 



2 ' 

2 arcsinh 



2 5 *\ 37 2 I Si} 

for x > | + s, 



(16) 



sinh ; 



cn(rea;cosh ^^,1/ cosh 4=p-) 



for < x < | - s, 



2 arcsinh 



sinh 






en 


e a ft(|a;|— x{) cosh 


H 


,1/cosh 







2a, for | - s < x < % + s, 



(17) 



4 arctanh |exp — n (j 
where we introduced 



tanh%} 



for x > j + s, 



A( 



- y 2 cosh (f) — 2 cosh i 



2 sinh —i/l 
2 V 



sinh — — / sinh — 
2 1 2 



A, 



(18) 

/ 2cosh(^ + 2a) -2cosh(^i + 2a), (19) 

cn(tt, k) is the Jacobi cosine-amplitude elliptic function with modulus £;[22, 23|, x\ is the inversion 
point where the electric field vanishes, (j)'{x\) = 0, and the electrostatic potential <p\ = <p(xi) is an 
integration constant to be determined by the boundary conditions (13) and (14). 

Matching the electrostatic potential 4>(x) at the planes x = ^ ± s by imposing the boundary 
conditions (13) gives 



2 arcsinh 



sinh 



cn 



K(|-s)cosh%i,l/cosh^ 



KS 



2e" Q cosh + a) 



■>XD.1 H J < 



cosh Sjs 



+ 



cosh 



(*+«)' 



where we introduced 

!F(<f) m ,(f)) = F ^arccos 
^• a (0i, 4>) = + 2a, + 2a) 



sinh — / sinh — 
2 7 2 



, 1/ cosh -y 



F |arccos 



sinh ( — + a ) / sinh ( — ha 



, 1/ cosh ( — + a 



(20) 
(21) 
(22) 

(23) 
(24) 



and .F(V') k) — J d0/ vl — /c 2 sin 2 # is the elliptic integral of the first kindf^, 23]. 

o 

On the other hand, matching the electric field <f>'(x) at the planes x = \ ± s by imposing the 
boundary conditions fll4| ) leads to 



e a [cosh(0< + 2a) — cosh(0> + 2a)] = cosh</>< — cosh< 



2 sinh 2 

2 



cosh(0i + 2a) = cosh(<£> + 2a) - 2e 2 ° sinh 2 



(25) 
(26) 
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Equations (21) an d (|25| ) represent a pair of coupled equations which can be solved for cj) m and 0>, 
since we can use ([20D and (p6|) to eliminate 4> < and (pi, respectively. Once obtained 4> m and 0> 
which solve equations and (p^), the electrostatic potential 4>(x) can be obtained by replacing 
them into the closed expression (\l7\j. To illustrate typical profiles for the non-overlapping regime, 
in Figure (|l]) we show the reduced electrostatic potential (j)(x) and the density profiles n±(x) for a 
fixed value of e, s and h. We also present the particle-density excess over the reservoir, 



and the charge density, 



n(x) = n + (x) + n_(x) — 2nn, 



p(x) = n + (x) — n^(x). 



(27) 



(28) 



The total free-energy density associated with the electrostatic potential (17) and the microion 
profiles (|8|) and (|9|) is obtained by replacing their closed forms into the total free-energy functional, 
given by (||), and performing the integrations. After some algebra, we obtain 



— / = 8ks ( 1 - 

-8^ r ((/> m ,(/) < 



2a 



2A( 



'in, V<) 



i< - 4coth — ) + 165 (<rt. M - - , 



cosh 

2 



sinh 2 ^sa 

^ + 2e- a A Q (0 i ,0<; 

cosh 



4coth 



+ a 



-2e~ Q A Q ,(0 i , 4>>) 



4coth 



+ a 



+ 16e" a [£, 



)+£«( 



)] cosh 



+ a 



-8e" Q [.F £ 



)+F a (<k, <!>>)] 



sinh' 



cosh 



+ a 



+ 4 sinh — 



4tanh — ^ 



where we introduced 

£($> m ,(/><) = £7 ^arccos 
£ a {<t)iA) = £{(t>i + 2a,(j) + 2a) 



sinh —51 / sinh — 
2 7 2 



, 1/ cosh — ) , 



E { arccos 



sinh ( — + a ) / sinh ( — ha 



,1/cosh ( -j 



a 



(29) 



(30) 



(31) 



and E(if>,k) = jd0Vl-k 2 sin 2 # is the elliptic integral of the second kind]22|, |23| . The closed 
o 

analytical expression (^9|) was checked against numerical integration of the free-energy density (H) . 
At the end of Section 4, in Figure (||), we present the total free-energy density as a function of the 
separation of the surfaces h for a fixed value of the depletion strength e and several values of the 
depletion range s. 

The disjoining pressure II is given by the negative derivative of the total free energy with respect 
to the separation of the surfaces, h, for constant depletion strength, e, and range, s, 



0U = 



"dh 



(32) 



where we introduced the dimensionless distances, h = nh and s = ks. After a lengthy calculation 
(see Appendix), we obtain a very simple final expression for the disjoining pressure, 



f3U = 4n sinh 2 



n [x 



0). 



(33) 
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The above simple analytical expression was checked against numerical differentiation of the free- 
energy density for the non-overlapping regime (|29|). Thus, it turns out that the disjoining pressure 
for the non-overlapping regime is given simply by the excess osmotic pressure of the microions at 
the midplane over the bulk (reservoir) pressure. Although it might be tempting to attribute this 
simple result to the contact-value theorem for charged plates p3, 25, E6|, we stress that this is not 



the case. Actually, an expression for the particle-density excess similar to the charged-plates case, 

n (x)=n(x = 0) + ^ [<//(£)] 2 , (34) 



holds only for < \x\ < 



cations n+(x) upon crossing the surfaces at x = =t ^ ± s, 



. Since there are non-vanishing discontinuities for the density of 

h 



An H 



h 

s 

2 

h 

2 +S 



= / /-|_ ^ x | 

'+• i x T 



h 

s 

2 

h 

2 +S 



n + [ x I 

x i 



n 



h 

s 

2 

h 

2 +S 



(35) 
(36) 



the corrected expressions for the particle-density excess for \x\ > 



arc 



n(x) 





h 








) 




2~ S 





n(x = Xi ) + ^ [<P'{x)] 2 , 



for 



h 




h 


s 


< x < 


2 +S 


2 







n(x) = n(x = 0) + An+ I x 





h 


^ + An+ (x = 


h 




^x = 




2 +S 


) 




2~ S 







n 



for x > 



+ s 



(37) 



(38) 



where we used (25) and (^) to simplify the above expressions. However, these additional con- 
tributions cancel when we evaluate the disjoining pressure (for the non-overlapping regime), and 
we obtain d33|), a result similar to the contact-value theorem expression for charged plates. Since 
the disjoining pressure ( p3| ) is always positive, the interaction between the surfaces for the non- 
overlapping regime is always repulsive. At the end of Section 4, in Figure (|j), we present the 
disjoining pressure as a function of the separation of the surfaces h for a fixed value of the deple- 
tion strength e and several values of the depletion range s. 



4 The overlapping regime, h < 2s 

In the overlapping regime, which occurs when the separation between the surfaces is smaller than 
the range of the depletion forces, h < 2s, the depletion zones associated with the two interfaces do 
overlap, and the generalized PB equation reads 



& 2 <j>(x) 
dx 2 



e 4a K 2 sinh [<f>(x) + 4a] , 



for < x < s - §, 



e k sinh [4>{x) + 2a] , for s 



(39) 



k 2 sinh ( 



for x > s + |. 
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The calculation is analogous to the case when there is no overlapping of the depletion zones, 



h > 2s. Now the pair of coupled equations to be solved for 
is given by 



us 



^2a 



+ 



= fax = 0) and </>> 
Fa(<t>U &<) 



(X = s 



e~ 2a cosh (if + 2a) 2e~ a cosh + a 



(40) 



and 



e 2a [cosh(0< + 2a) - cosh(0> + 2a)] = e [cosh(c/>< + 4a) - cosh(0 m + 4a)] - 2 sinh 
where (f) < = fax = s — ft) and fa = fax = x{) are eliminated by using the relations 



2 v> 



(41) 



2 arcsinh 





sinh 


O + 2«) 




cn 


e-** (s - f ) 


cosh 


[if + 2a) 


,1/cosh (if + 2a) 



cosh^i + 2a) = cosh(0> + 2a) - 2e 2Q sinh 2 



4a, (42) 
(43) 



Once solved the system ([40|) and (pl|), the electric field and the electrostatic potential can be 
obtained by replacing the solution (c/> m ,^>) into the closed expressions, 



Ke 2a A 2a [4> m ,fax)} 



for < x < s - §, 



<f>'(x) = < Ksign (x - xi) e a A a [fa, fax)] , fors-|<x<s + | 



2 ; 



(44) 



-2k sinh 



2 ' 



for x > s + ft, 



0(x) 



2 arcsinh i 
2 arcsinh 



sinh(%i+2a) 



cn c ™ri 



cosh ( if +2a) , 1 / cosh ( if +2a) ] 
sinh 



4a, 



for < x < s - §, 



e a ft(|a;|— x\) coshl -f+ct ) ,1/coshl -f+ct 



2a, for s-f <x<s + |, (45) 



4arctanh |exp 



-K IX 



tanh^} 



for x > s + 



with the inversion point given by 



^2a 



F a (fa 



, H J <, 



-2a 



cosh (if 



+ 2a e" a cosh 



(| + a) 



(46) 



To illustrate typical profiles for the overlapping regime, in Figure (||) we show the reduced elec- 
trostatic potential fax) and the density profiles n±(x),n(x) and p(x) for a fixed value of e,s and 
h. 

Again, it is possible to obtain a closed analytical expression for the total free-energy density, 



=M-i)( 



1 - e 



-Let 



+ 



4h( 



l-e~ 2a ) +— / aux , 
n 



(47) 



n 



a „ x = 2e 2Q A 2Q (</> m ,<A<) 



!>< -4coth ( -y + 2a 
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+ 1 <ie- 2a £ 2a (0m, <A<) cosh ( ^ + 2a ) - 8e~ 2a ^ 2l , ( 



sinh 2 (%- + 2a 
cosh + 2a 



-2e- Q A Q (</- i ,0 < ; 



; — 4 coth [ + n 



+ 2e _Q A 



' — 4 coth [ -I- n 



-16e a [£ a (<&, 0<) - </>>)] cosh f -j + a 

+8e- Q [^ a (</» i , ( /> < )-^( ( /» i ,0 > )] ' 2 



+ 4 sinh — ^ 4 tanh ^ 



) , (48) 



cosh(^ + a) 2 V r " " 4 

which leads, after some algebra (see Appendix), to a simple expression for the disjoining pressure, 



2n 



nix 



e 4a cosh( 



+ 4q) - 1 + 2e 



-<t>< — 6a 



sinh 2q 



0) + 2An+(x = s 



2' 



(49) 



The above simple analytical expression was checked against numerical differentiation of the free- 
energy density for the overlapping regime ([47|). It should be remarked that, in this case, the 
disjoining pressure has not the form of the expression given by the contact- value theorem for 
charged plates ^4|, |2(|. An additional contribution due to the discontinuity An + (x = s — |) of 
the density of cations upon crossing the surface located at x = s— |, appears. Contrary to the non- 
overlapping regime, this additional contribution does not cancel when we evaluate the disjoining 
pressure. According to this imbalanced pressure acting onto the neutral surfaces, this leads to an 
effective attraction between them. 

Figures @ and (Q) show the total free-energy density and the associated pressure as a function 
of the separation of the surfaces h for a fixed value of the depletion strength e and several values 
of the depletion range s. The nature of the interactions changes from attractive to repulsive at a 
separation h = 2s. 



5 Concluding remarks 

We have proposed a new mechanism for attraction between neutral plates immersed in a monovalent 
electrolyte solution, which does not include any correlation or fluctuation effects. The electrostatic 
potential and the density profiles of the microions are obtained from analytical solutions of the 
generalized PB equations, which include non-electrostatic depletion interactions. Explicit analytical 
expressions of all thermodynamical properties, including the disjoining pressure, were obtained. 

We found that the repulsive interactions at large separations become attractive when the sep- 
aration between the plates is decreased. The range of the attractive forces is closely related to the 
range of the non-electrostatic depletion interactions introduced in the formulation of the model. 
Although this result is not at all surprising, since the attraction is induced by the imbalanced 
pressure originated from the ionic depletion in the region between the two approaching surfaces, 
we found that the disjoining pressure has not the form of the expression given by the contact-value 
theorem for charged plates. 

The proposed mechanism could mimic neutral surfaces immersed in an electrolyte solution 
containing ions of different sizes. We expect to observe attraction when the separation between the 
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two surfaces is comparable to the size of the smaller ions. We stress the fact that we do not include 
any non-mean-field effects to obtain attractive forces. Surely, for short separations, other features 
should be taken into account, as for example, the discreteness of the charges and ionic correlations. 
However, using an exactly solvable model, we showed that the inclusion of non-mean-field effects 
is not a necessary condition to obtain attractive interactions. 



A The disjoining pressure 



The disjoining pressure II is given by the negative derivative of the total free energy with respect 
to the separation of the surfaces, h, for constant depletion strength and range, e and s, 



n 



no dh 



k df dfa 



no d(j) m dh 



k d/ d(f) < 
no d(j) < dh 



k df dfa 
no d4>i dh 



« df d4>> 
no d(j)y dh 



(50) 

where we introduced, for convenience, the dimensionless distances, h = nh and s = us. The (four) 
derivatives of the free energy which appears into (|50T) , df/d(p, with <p = ((f) m , 0<, fa, <j>>), can be 
obtained directly by using the free-energy expressions (29) and (^). On the other hand, the partial 
derivatives d<p/dh\ ai s are obtained by the matrix product 



d(p < dfa dfa 



dh dh ' dh' dh 



d(h, s, u, v) 

d((/>m, <f><, fa, <f>>) 



(1,0,0,0), 



(51) 



where u and v are the boundary conditions ( |14| ) written in a parametric form involving <p m , cj) < , fa 
and 0>. We will give the explicit expressions of u and v (and their derivatives) when we treat 
separately the non-overlapping and the overlapping regimes. 



A.l The non-overlapping regime 

For the non-overlapping regime, the total free-energy density is given by (|29|), with derivatives 



K df 



no dfa 



4£(^ m ,0<)sinh- 



2 sinh ( 



n dfa 
k df 2 sinh (j) < 



2 A( 



J m, <+>< 



2 tanh 2 coth \ , 
2 2 J 



^-2tanh^ + 2e " aSmh ^< + 2a ) 



A 



nod(/>< A((f) m ,(/) < 

— TT = 4e ~ Q ^<) + S <*(&> ^>)] sinh (t + a 

n dfa \ 2 

2e~ a sinh (fa + 2a) 



)< — 2 tanh ( — — \- a 



2e~ a sinh (fa + 2a) 
A Q ((/>i,0>) 
K df 2e~ a sinh (</>> + 2a) 



< — 2 tanh 2 [ — ^ + a ) coth ( ^ + a 



> — 2 tanh 2 ( + a I coth ( — + a 



i> — 2 tanh ( — — h a 



+ 2 cosh — [<j) > —2 tanh 



The defining equations for h,s,u and i> are 



^ _ 2F(4>m,4><) , ^a{fa,4><) +J r a(fa,4>>) 

cosh if e- a cosh(^ + a) 



(52) 
(53) 



(54) 
(55) 

(56) 
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_ Fgjfa, &<) + Fajfa, <!>>) 

2e-°cosh(^- + a) 
u = e~ a A a (<&, 0<) + A(0 m , </><) = 0, 
v = e- a A Q (0 i ,0 > ) + 2sinh^ = 0, 



(57) 

(58) 
(59) 



and their derivatives necessary for the evaluation of the Jacobian d(h, s, u, v)/d((p m , (p < , 4>u 4>>)i 



ds_ 

dji 

ds 

dj; 

du 
d(j) m 
du 



Si 



4m, 0<) 2 tanh 3f coth ^ 



sinh : 



+ 



1 



d/i 



d/j 



d</>< A(0 m ,^<) e- a A a (0i,0<) 

dfr _ [£g (<&,</><) +g Q -(0j,0 >)] 

d0i 2e- a sinh(^ + a 

d/i 1 



coth 



e" a A 



+ 



coth + a 
e- a A a (<k,(l>>) 



d<£> e-«A a (0i, </>>)' 
ds 



0, 



1 



ds 



d0< 2e- a A a (0i,0<) 



(60) 
(61) 



tanh ( y + a ) , (62) 



(63) 
(64) 
(65) 



4e" a smh +a) 
1 



coth (if + a) coth + a) 



2e- a A Q ( ( /> i 



2e- a A a (0i,0>; 



2e" Q A Q (0i, 0>) 
sinh (6m 



A( 



sinh 0< e a sinh (0< + 2a) 



+ 



du e _a sinh (0; + 2a) 
d0i Aq^i,^) 
du dt> di> 

i 7 



d^> d^ m d<^. 

di> e _a sinh (0; + 2a) 

dcf)i A a ((f>i, (f>>) 

dv e~ a sinh (</>> + 2a) 



+ cosh ■ 



tanh ( — + a ) , (66) 



(67) 
(68) 
(69) 
(70) 
(71) 
(72) 
(73) 



d<£> Aa,(0i,0>) ^ 

Putting all together into the expression for the disjoining pressure (|50|), leads to a very simple 
final result, 



/5II = 4?t-o sinh 2 = n(x = 0). 



(74) 



Although it might be tempting to attribute this simple final result to the contact-value theorem 
for charged plates, this is not the case (see discussion at the end of Section 3). 
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A. 2 The overlapping regime 

For the overlapping regime, it is convenient to apply the parametric differentiation just on the last 
term, / aux , of the total free energy (|47|), since the two first terms yield a constant contribution to 
the disjoining pressure, 



n 



no dh 



2 (l- e ~ 4Q ) -4(l-e" 2Q ) 



« dfa. 



n d</> E 



n &<t>< 



l<--°£ 2 a(&n,</k)sinh( Z f + 2a 



The derivatives of the last term 

2e~ 2a sinh (0 m + 4a) 

A2a(0m,0<) 

2e" 2a sinh (0< + 4a) 

A 2 a(</>m,</><) 

2e" a sinh (</>< + 2a) 



no <9/i 

of the total free energy, / aux , are given by 



2 tanh 2 — + 2a coth — + 2a 



2tanh 



+ 2a 



< — 2 tanh ( — + a 



K d/ aU5j 
no dcpi 



+ 



Ao,(0i,0<) 

4e" Q [5« (0i, 0<) - </>>)] sinh 

2e _a sinh + 2a) 



A a (^i,^><) 
2e" a sinh (<fc + 2a) 



2 tanh^ 



+ a 



— + a ) coth I — ^ + a 



W d/ a 

no d(j> 



2e~ a sinh (0> + 2a) 



> — 2 tanh 2 ( — + a ) coth ( — + a 



> — 2 tanh ( — — h a 



A a (0i,0>) 

Now the defining equations for h,s,u and w are 

F a {(t>h 4>>) ~ F<x{4>h <t><) 



+ 2 cosh — [(/>>— 2 tanh — 



e~ a cosh ( 2 + Q 



-F 2 



+ 



e~ 2a cosh (%i + 2a) 2e" a cosh + a) 



u = e Q A Q Oi,i 



-2a 



A 



2<H<Pm,<P< 



) = o, 



v = e- a A Q (<Ai, </>>) + 2 sinh ^ = 0, 



(75) 



(76) 



(77) 



(78) 
(79) 

(80) 

(81) 

(82) 
(83) 



and their associated derivatives for the evaluation of the Jacobian d(h, s, u, v)/d(cf> m , (/><, 0i, <f>>), 
dh 



d(p m 
dh 



0. 



1 



e" ft A 



(84) 
(85) 
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Figure 1: Reduced electrostatic potential (j)(x) and density profiles as a function of the distance 
x for the set of parameters ne = 1, us = 1/2 and nh = 5 (non-overlapping regime). All densities 
are normalized to the salt reservoir density, uq. The positive portion of the particle-density excess 
n(x) was amplified by a factor of 25. 
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Figure 2: Reduced electrostatic potential 4>(x) and density profiles as a function of the distance x 
for the set of parameters ne = 1, ks = 1/2 and nk = 1/2 (overlapping regime). All number densities 
are normalized to the salt reservoir density, uq. The positive portion of the particle-density excess 
n{x) was amplified by a factor of 10. 
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0.5 1 1 — — -> 

2 4 6 

Figure 3: Reduced total free-energy density / as a function of the separation of the surfaces 
h for a fixed value of the depletion strength (kc = 1) and three values of the depletion range 
(ks = 1/4, 1/2, 1). Although the free-energy density itself is continuous upon crossing the separation 
h = 2s, it has a kink at this special value, giving rise to a change between attractive and repulsive 
forces (see next graph). 
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Figure 4: Disjoining pressure II as a function of the separation of the surfaces h for a fixed value 
of the depletion strength (ne = 1) and three values of the depletion range (ks = 1/4,1/2,1). 
This graph corresponds to the negative derivative of the curves of the previous figure. Note the 
discontinuity of the pressure at the separation h = 2s, associated with the kink of the free-energy. 
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